AD-A070  215 


UNCLASSIFIED 


WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/G  12/1 

A COMPARISON  OF  SAMPLIN6  PLANS  FOR  BAYESIAN  ESTIMATION  OF  A POI— ETC(U) 
MAY  79  R L WARDROP  DAA629-75-C-0024 

MRC-TSR-1956  NL 


I I 

*§070215 


MRC  Technical  Summary  Report  #1956 


A COMPARISON  OF  SAMPLING  PLANS 
FOR  BAYESIAN  ESTIMATION  OF  A 
POISSON  PROCESS  RATE 


Robert  L.  Wardrop 


UNIVERSITY  OP  WISCONSIN  - MADISON 
MATHEMATICS  RESEARCH  CENTER 


r" 

L 


A COMPARISON  OF  SAMPLING  PLANS  FOR 
BAYESIAN  ESTIMATION  OF  A 
POISSON  PROCESS  RATE 

Robert  L.  Wardrop 
Technical  Summary  Report  #1956 
May  1979 
ABSTRACT 

Bayes  estimation  of  the  arrival  rate  of  a Poisson  process  is  studied  in 

A 2 — p 

this  paper.  For  any  loss  function  in  the  family  L^  = (0-0)  0 , < p < «•, 

a simple  sequential  procedure  is  introduced  which,  based  on  the  criterion 

of  minimizing  expected  cost  (estimation  error  plus  sampling  cost) , is  either 
optimal  or  asymptotically  optimal.  The  procedure  is  compared  to  Type  I 

and  II  censoring  - the  comparison  should  be  useful  to  experimenters  choosing 
between  the  three  sampling  plans. 
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SIGNIFICANCE  AND  EXPLANATION 


When  trying  to  explain  or  analyze  events  that  occur  randomly  in  time  or  in 
space,  one  tends  first  to  test  whether  the  events  are  governed  by  a Poisson  dis- 
tribution. The  classic  example  concerns  the  number  of  soldiers  killed  per  year 
from  the  kick  of  a mule  in  the  Prussian  Army  in  the  early  1800's,  and  applications 
have  continued  to  this  day  in  many  different  contexts,  military  and  otherwise. 

Usually,  the  mean  arrival  or  occurrence  rate,  0,  of  a Poisson  process  is  un- 
known. This  paper  derives  optimal  and  approximately  optimal  procedures  for 
sampling  from  a Poisson  process  and  estimating  the  value  of  0.  Three  classes  of 
procedures  are  considered.  First  is  "Type  I censoring"  in  which  the  length  of  time 
the  Poisson  process  is  observed  is  fixed  in  advance.  Second  is  "Type  II  censoring" 
in  which  the  number  of  arrivals  or  occurrences  observed  is  fixed  in  advance.  The 
final  class  of  procedures  is  "sequential"  in  the  sense  that  neither  the  lenqth  of 
observation  nor  the  number  of  occurrences  observed  are  fixed  in  advance.  Instead, 
the  outcome  of  the  process  prior  to  any  point  in  time  may  be  used  to  decide 
whether  to  continue  observing  the  process  beyond  that  time'. 

The  best  procedures  in  each  of  the  three  classes  are  compared.  The  results 
are  useful  to  the  experimenter  who  wants  to  choose  an  efficient  sampling  plan. 

Accession  For” 

NTIS  U.iAAI 

DDC  TAB 

Unannounced 

Justification 

By 

Dlst  r l1"  t.  t cn/ 

Avail r>  h J ’ 1 1 y_  C o das  I 

Aval  l and/or 

Dlst  special 

[A  I 

Tit  t osponsibi 1 ity  for  the  wording  and  views  expressed  in  this  descriptive  summary 
with  MRC,  and  not  with  the  author  of  this  report . 


A COMPARISON  OF  SAMPLING  PLANS  FOR  BAYESIAN  ESTIMATION 
OF  A POISSON  PROCESS  RATE 


Robert  L.  Wardrop 


Section  l:  Introduction  and  Notations 

Conditional  on  the  value  of  ft  > 0,  let  X(t)  , t j>  0,  be  a Poisson  process  with  arrival 

rate  0.  Sot  t()  « 0,  and  for  i » 1,2 let  tj  « inf{t  ! X(t)  - i ) be  the  time  of  the 

ith  arrival.  Bayes  estimation  of  0 will  be  studied  in  this  paper. 

Assume  that  the  loss  incurred  from  estimating  ft  by  6 is  qiven  by 


(1.1) 


L^lft.ft)  - (ft-ft)2e”P 


for  some  p,  -«  < p v •».  loss  functions  of  the  form  (1.1)  have  been  proposed  by  numerous 
authors,  including  Hodqes  and  Lehmann  (19S3)  (p  « 1),  Ovoretzky,  Kiefer  and  Wolfowitz  (19>1) 

(p  - 1),  El-Sayyad  and  Freeman  (1973)  (p  - 0,1,  or  2),  Novic  (1977)  (0  ^ p £ 3)  , and  Shapi>o 
and  Wardrop  (1977,  197H)  (0  p <_  3).  These  papers  present  some  justification  of  various 
choices  of  p,  especially  p • 0,1,  or  2.  Results  will  be  obtained  for  all  real  p because 
this  does  not  increase  the  difficulty  of  the  proof s and,  more  importantly,  because  it  provide- 
additional  insiqht  into  the  behavior  of  the  sampling  rules  considered. 

Assume  0 has  prior  distribution 


. a ft,-,-! 

X (0)  - I'("0>  % ft  exp(-B08) 


with  t<()  ' 0 and  a • 0 v p<.  (Many  of  the  results  obtained  are  true  with  a «.  p,  thi  v>  i 1 1 


be  discussed  uqain  in  Sections  3 and  4.)  Denote  this  distribution  by  ftn^,!^)  . For  ( 


let  .*’(t)  denote  the  s iqma-a  Iqebra  of  events  generated  by  IX(s),  0 - s • t).  The  postoii 


distribution  of  ft  qiven  F(t)  is  f(a(t)  ,ft(t) ) , «(t)  - o 1 t X(t)  and  ft(t)  * » t. 


loss  l.  , the  Bayes  estimator  of  ft  given  V(t)  is 


(1.2) 

and  the  |vintei  ioi  i’*|»it ed  loss  Is 


ft)((t)  - (a(t)-p)  (ft(t) ) 


-l 
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(1.3)  E(L  (0,6  (t))|F(t,(  = 6(t)P_2r(ot(t)  + 1 - p)  r (a  ( t) ) -1  . 

p p 

It  can  be  shown  that  (1.2)  and  (1.3)  remain  true  with  t replaced  by  o,  a stopping  time  with 

respect  to  F(t)  , t ^ 0,  if  F(o)  is  qiven  its  usual  definition  (Shapiro  and  Wardrop  (1977)). 

For  a precise  method  to  discourage  samplinq  indefinitely,  let  c > 0 be  the  cost  of 

A — 

observing  one  arrival  of  the  process,  let  c^,  0 be  the  cost  of  observing  the  process  for  one 

unit  of  time,  and  assume  that  cft  v c^  > 0.  The  total  cost  of  observing  the  process  for  t 
units  of  time  is  defined  as 

(1.4)  C(  t)  = C (t)  « 6 (t)  P_2r  (a(t)  + 1 - plTfaft))'1  + c X(t)  + ct  . 

p AT 

Different  sampling  plans  will  be  compared  on  the  basis  of  expected  total  cost. 

To  motivate  later  results,  note  that  the  values  of  c and  c represent  the  cost  of 

A T 

sampling  measured  in  the  units  of  the  loss  function  L (this  is  clear  from  the  definition  of 

P 

C) . Intuitively,  asymptotic  results  should  be  obtained  by  letting  the  cost  of  sampling  relative 
to  the  cost  of  estimation  error  decrease, because  this  will  encouraqe  longer  observation  of  the 
process.  One  way  to  achieve  this  is  to  let  and  c tend  to  0.  Another  approach  (follow- 

ing El-Sayyad  and  Freeman  (1973))  would  be  to  define  total  cost  as  Oft)  « DE(Lp|F(t))  + 
cftX(t)  + c t and  let  D ■+  °°  while  holding  cft  and  cT  fixed.  The  two  methods  are  obviously 
mathematically  equivalent;  in  this  paper  the  first  method  will  be  used. 

The  total  cost  function  may  be  written  in  two  other  ways  which  will  be  useful  later: 

(1.5)  C( t)  = E(01_p|F(t))6(t)_1  + cftX(t)  + c t , 


(1.6)  C(t)  = E ( 02~P | F( t) ) (a( t)  + 1 - p)'1  + cftX(t)  + c^t  . 

For  t > 0,  let  Y(t)  = EfB1  P|Ft>  and  Z(t)  = E(02  P|F(t)).  By  a well  known  theorem 
{ Y( t)  , t 2,  0}  and  (Z(t),  t 0}  are  uniformly  integrable  martingales.  Thus  Y(«>)  * lim  Y(t) 

t-K® 

and  Z(")  = lim  Z(t)  exist  almost  surely  and  Y(»>)  = 0 P and  Z ( <*>)  =0  P.  For  b > 

^ -MO 


(1.7) 


vb  ■ E‘°b)  * r<vb)r(v'1C  • 

In  Section  2 a stopping  time  t will  be  defined  and  shown  to  be  either  opt  w)  or  asymp- 
totically optimal  for  all  p,  cft,  c^,,  and  6^.  In  addition,  the  limitinq  form  of  E(C(t)) 

will  be  given. 

In  Section  3, nonsequential  sampling  plans  will  be  considered,  namely  Type  I censoring 
(observing  the  process  for  a predetermined  length  of  time)  and  Type  II  censoring  (observing  the 
process  until  a predetermined  number  of  arrivals  are  observed).  Using  the  criterion  of  mini- 
mizing expected  total  cost,  the  Bayes  Type  1 censoring  (Bl)  and  Type  II  censoring  (B2)  proced- 
ures are  obtained  explicitly  along  with  their  respective  expected  costs  (V^  and  V2> . The 
values  Vj  and  V.,  are  compared  asymptotically  to  determine  the  cases  in  which  Bl  is 
superior  (inferior)  to  B2. 

In  Section  4,  the  results  of  Sections  2 and  3 are  combined  to  determine  how  much  better 
t performs  than  Bl  or  B2.  An  explicit  asymptotic  measure  of  the  improvement  is  given. 


Section  2t  Sequent ial  Sampling  Plana 

All  stoppinq  times  are  with  respect  to  (F(t),  t > Ol.  The  stoppinq  time  o is  called 
optimal  if,  and  only  if  E(C'(o))  * inf  E('’(r))  with  the  infimum  taken  over  all  stoppinq  times 
0.  In  such  problems,  it  is  often  useful  to  compute  the  infinitesimal  generator  of  the 
stochastic  process  C(t).  It  is  defined  to  be 

*'(t)  - lim  h"1|E(C’(t»h)  - f(t>|F(tin  • 

h+0 

Using  (1.5),  it  is  easy  to  show  that 

Ai'(t)  - -6(t) --2Y (t>  ♦ c altieitl’1  ♦ c , 

A T 

for  Y(t)  defined  in  Section  1.  Intuitively,  as  lonq  as  Af(t)  < 0,  sampling  should  continue 
since  the  total  cost  is  "expected"  to  decrease.  A natural  stoppinq  t lme  to  consider  is 

(2.1)  T « least  t > 0 such  that  AC(t)  > 0 , 


Y(t)(Mt)'2  < c.aftllMt)'1  ♦ c . 

- A T 

For  p an  inteqer,  the  rule  t is  easy  to  use;  for  example  if  p - 0,  t stops  the  first  time 

(2.2)  a(t)£(t)  2 < c o(t)  ♦ c_p(t) 

— A T 

The  left  side  of  (2.2)  is  the  posterior  variance  of  P while  the  riqht  side  is  approximately 
the  total  cost  of  sampling.  In  fact,  for  any  p,  t stops  the  first  time 

E ( L ( 9 , 9 ) ( F ( t ) ) cAl'<t)  + cT*(t)  , 

which  generalizes  the  above. 

The  following  result  on  the  optimality  of  t was  obtained  independently,  using  different 
meth.vls  of  proof , by  Novic  (1>»77)  and  Shapiro  and  Wardrop  ( 1 *»77 > . 


Theorem  2.1.  In  the  cases 

( i > 0 p < 1 and  c • 0 , 

_ T 

(ii)  1 < p < 2 and  all  c.,  c_,  or 

A T 


-4- 


(ill) 


1 ' 1 and  o » 0, 
A 


\ given  by  (2.1)  i * opt  imal  foi  all  ' 0 and  ^ p 


Foi  situations  not  oovcred  by  Theorem  . 1 , including  t lie  interest  lnq  case  p - 0 and 

e ' 0,  t he  optimal  stopping  t \me  is  not  Known . Moreover,  thorn  are  not  any  gennial  existing 

results  on  thr  limiting  form  of  tho  expected  total  cost  of  tho  optimal  rule.  In  this  scot  ion 

it  will  bo  shown  that  i is  asymptotically  optimal  for  all  p,  c , c . tt  and  a,('  p)  , and 

A T O 0 

tho  limiting  foim  ot  K(('(t)l  will  lw>  obtained.  First  a lower  bound  foi  tho  asymptotic  expect 

oil  total  cost  of  any  sequence  of  stopping  times  is  obtained  in  Theot  em 

l.emma  2.1.  For  l>,  V ' 0 random  variables,  min  F.(UV  ' ♦ V!  - 2K(t''  ’! 

V 

Proof ; For  x,  y ' 0,  g(x)  - yx  ' ♦ x achieves  a unique  minimum  of  .'y ' * at  x - v ' ", 
Therefore  F.FU’V'1  ♦ v|l'l  ' 2K(U1  ‘l  . 

Them  em  2.2. 


(l)  11  c * 0,  let  c - o(c  ) be  any  family  of  stopping  times,  then 


Iim  mf  o*'  ’ F ( > ' < c 1 ) ' 2KY  (<*> 1 * - 21‘ia,  ♦ (1-p)  2)  t»  *V'~  * ’ *V  iv»0>  ” * 

V° 


(i»)  It  ^ 0,  wiito  c*^  • * 0)  . l.oi  o - iMo  ) bo  any  family  of  stop\mui  » • m 


Iim  »nf  i*  ' N 2RiaYi*)  ♦ > * * - * ♦ 0*  *') 

» A — 


Vl  oof  : IVi  ( i ) , 


• 1/2  • t 2 - 1 \ 2 

\ *m  inf  <*  K “ 1 im  inf  K(Y(oK'  ♦ v'  K ' 

T T v'  T 

c_  ‘0  c >0 

t r 


2_  l \m  inf  Jt?Y(ot 
v'  *0 


U. 


p 


by  Lemma  2.1  and  the  fact  that  (Y(t)'  * , 0 *•  t < •*•)  is  a supermartinqale. 


For  ( x a ) 


-1  '2 

1 im  inf  E(  (o) ) 

CA*°’VaCA 


lim  inf 


v-  i (a8(o)  ♦a(o)  + l-p)  r (g(o)  tl-p) 
\ 1/2 
\ A 


(a6(o)+n(o)+l-p)8(o)2*pr(a(o)) 


1/2 

+ c,  (a(o) 

A 


♦ 1-p  ♦ a8(ci)  )\  * 


lim  inf  EUlv'1  + V) 


1/2 

with  U « aY(o)  ♦ Z(o)  and  V » (a8(o)  + a(o)  + 1 - p) . The  result  now  follows  as  in  (i) 


To  obtain  the  limiting  value  of  E(C(l)),  the  followinq  result  is  needed. 
Lemma  2.2. 

(i)  If  c » 0,  then 

A 

lim  ci/2ET  - E0(1-P)/2  . 

o*0  T 
T 

(ii)  If  c_  > 0,  and  c.  * ac.(a>0),  then 
A T A — 

...  . . 1/2 

lim  c ~E(a8(i)  ♦ a(il)  » E(aO  P ♦ 0‘  *') 

, A 
c.  *0 
A 

C -ac 
T A 

2 

Proof : In  case  (i),  from  the  definition  of  t,  Y ( T > <_  0^8(1)  . Thus, 


(2.1) 


cV’EISIt))  ' E(Y(l)  1 ') 
T - 


For  » - 0,  on  the  set  T - • , Y(t-r)  ' ) * » yielding 


L_ 


V'T  Etjltt-.))  v E(Y(l-.)  ) . 

tt  is  oj»y  to  sm  that  a*  o * 0,  i » • and  Y(t)  * Y(*l  with  probability  on«.  Moreover, 

both  Y(t-*)1  * and  Y < 1 1 ' * are  bounded  above  by  (sup  Y ( t ) > ' *.  This  latter  random  variable 

t>0 

is  integrable  because 


P ( ( sup  Y ( t ) ) ' a)  - lis  p(  sup  Y(t)  > a ) s a B(Y  (> 


r *-'■  o<t<T 


Part  111  now  follows  fieai  l-'.'l  and  1 ^ . -4  > . 


for  ease  (li),  the  definition  of  l gives 


aY(t)  ♦ 2 (i)  e^(atl(  r)  ♦ a(i))*  » o ( 1-p)  (ati(  t ) ♦ a(t)) 
and  on  the  set  t » , > 0, 

2 

aY(i-t)  ♦ T.(i-r)  - e (att(t-c)  ♦ a(t-t))*  ■*  c.  (1-p)  (att  ( t -« ) + o(t-t)) 

A A 

The  lemainder  of  the  proof  is  similar  to  case  (i)  and  will  not  be  given. 

The o^l  em  . 1 , 

(i)  tf  c * 0,  then 


lim  c *E  ((*(  t ) ) * -K  ( 0 1 1 l'  '1 

c »0 

T 


(ill  tf  c ' 0,  and  » ae^ianll  , then 


. , . , 1/2 
e,  *K(t'(tl)  - 2E(aO  l'  ♦ tl*'Pl 


Pi  tH't  : Foi  ( i > , 


M\d  toi  ( \ i)  , 


l'{"  ■ V,  4 V --  -v.  * 


.'(r>  • Y(i)tMt)  ♦ t'  (At  ♦ X ( i ) ) ' (atMO  ♦ 

A A 


-7 


The  ile  si  red  results  follow  from  Theorem  2.2  and  Lemma  2.2. 


In  view  of  the  results  presented  in  Theorem  2.2  and  2.3,  say  that  \ is  asymptotically 


-1/2 

optimal.  Note  that  if  both  e and  c are  positive  then  1 lm  c F ( % ( t ) ) must  be  com- 

A * .A 

c *0 

A 


pu  t ed  nume  i i c a 1 1 y . 


1 


Section  \ . 


Types l and  II  Censoring 

By  Theorem  2.1,  the  optimal  sequential  procedure  is  type  I censoring  if  either  * p 

or  c -»  0 and  p 1,  and  it  is  type  IT  censor inq  if  either  c * 0 and  p *2  or  c 
A t A 

and  p * 3 (the  cases  with  p * 1,2  were  obtained  by  El-Sayyad  and  Freeman  097})).  in  some 
applications  it  may  not  be  feasible  to  use  r,  and  the  experimenter  must  choose  either  Type  1 
or  II  censoring.  The  results  of  this  section  should  be  helpful  in  making  that  choice. 

The  Bayes  type  I censoring  procedure  (Bl)  is  that  t*  which  minimizes  E(J(t))  over  all 
t 2,  0.  Using  representation  (1.5)  for  C( t)  , it  is  easy  to  verify 

1/2 


t*  - {(v  (c  v ♦ c : 

l-p  A l T 


-1, 


V 


and  for  t*  > 0, 


o.l)  Vj  - E(C(t*))  - 2lVP<Vl  + CT,]  ' VVl  + Ct’  ' 

with  v»  given  by  (1.7). 
b 

The  Bayes  type  II  censoring  procedure  is  that  integer  n*  which  minimizes  E(C(t^))  ov«  » 
n * 0,1,2,...  . Treating  n as  a continuous  variable  and  using  representation  (1.6)  foi 

/(t),  it  is  easy  to  show 

1/2 

U.2)  n*  » { lv2_p <ca  + CTV-1)  ' * (a0  + 1~)') '* 

providevi  that  either  a_  > 1 or  c_*  0.  If  a_  < 1 and  c s 0,  then  n*  = 0 since 
0 T 0 — T 

c F(t.)  Foi  the  remainder  of  this  section  assume  that  either  a **  1 or  c =0.  It 

T 1 0 T 

n*  > 0,  then  set 


(*•3)  V F.(i’(t  ))  * 2[v  (c  + c v . )1*  " - (a  ♦ l-p)  (c.  + c v ,) 

2 n*  2-p  A T-l  0 r A T - 1 

For  * he  remainder  of  the  section  assume  c and  c are  small  enough  to  insure  n* , 

A T 

t*  > 0. 

Remark  * . 1 . The  value  n*  given  by  (3.2)  is  not  necessarily  an  integer,  so,  strictly  speaking, 
B2  is  either  ln*|  oi  [n*  ♦ 11,  1*1  the  greatest  integer  function.  Also,  the  expected  ‘otal 


cost  of  B2 


is  not  exactly  V2  if  n*  is  not  an  integer.  For  qiven  values  of  c»0,  BQ,  p, 
cft  and  cT,  one  may  compare  and  (or  V2 ' s 'exact'  version)  as  an  aid  in  choosing 

between  B1  and  B2.  In  the  remainder  of  the  section  and  V2  will  be  compared  asymptotic- 
ally as  sampling  costs  go  to  zero  (the  effect  of  n*  not  being  an  integer  disappears  in  the 
limits  . 

Comparison  of  and  V 2 when  cft  = 0. 

If  p = 2,  then  V,  = V.  for  all  c >0.  Define 
12  T 


(3.4) 


Vp'"0>  = 


lim 

c -K> 
T 

cA=° 


V V 
12 


= Iv, 


-1  -1  . 

V V 

-p  -1  2-pJ 


1/2 


= ( (ciy-1)  (aQ+l-p)  ] 1,  for  aQ  > p v 1 , 

by  (3.1),  (3.3)  and  (1.7).  It  is  not  difficult  to  show  that  result  (3.4)  remains  true  for 

a > (p-1)  v l.  Moreover,  if  0 v (p-1)  < a < 1,  then  R. (p,a„)  = 0.  Clearly  R. (p,a  ) < 1 

U A 0 n U 

if  and  only  if  p < 2.  In  words,  asymptotically  B1  gives  a lower  expected  cost  than 

B2  iff  p < 2.  As  shown  above,  the  asymptotic  expected  savings  in  using  B1  instead  of  b2 

can  be  100%  (when  R = 0) . Conversely,  the  asymptotic  expected  savings  in  using  B2  instead  of 

_!  !/2 

B1  can  approach  100%  when  R^  ■*  °°  (e.g.  take  p = 3,  otg  = 2 + e . e > 0;  then  Rfl=  [(l+e)e  ] -»  °° 
as  c ->-0).  Thus,  depending  on  the  values  of  p and  aQ  there  can  be  a tremendous  difference  be- 
tween and  V2 . If  (p , aQ)  is  near  unity  then  the  difference  (based  on  expected  total  cost) 

between  B1  and  B2  is  not  dramatic  and  the  experimenter  may  choose  to  use  the  procedure  which  is 
easier  to  implement. 

Comparison  of  V,  and  V,  when  c =0. 

1 2 T 

If  p = 1,  then  V,  = V,  for  all  c >0.  Define 
12  A 


(3.5) 


VP'V  = lim  V1V2X  = <Vl-pVlV2-p) 

c ->0  ^ ^ 

A 

c =0 
T 


1/2 


, -1 


1/2 


= (a0(aQ+l-p)  I , for  aQ  > p v 0 , 
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hy  (4.1,  (l.l)  and  (i.7).  it  is  not  difficult  to  show  that  (l.S)  remains  true  for 

c»0  ">  (p-i)  v o (for  B2  one  need  not  require  s 1 since  * 0).  Clearly  R^(p,a^)  1 

if  and  only  if  p < l and  R takes  on  all  values  in  (0,***)  (e.g.  R^,  ♦ * as  1 (p-l)  • 0i 

for  p * 0*  H,j,  0 as  n -*  0).  The  discussion  qiven  alcove  for  » 0 is  also  televent 

in  this  case. 

Comparison  of  V,  and  V*  when  c -c  > 0. 

1 . A T 

Pot  a s 0 define 


( l.b) 


K<p'VaV 


lim  V V * 

c >0  12 

A 

c -ac 
T A 


<VatV  {(Vl-p,<Ua8o<V1,  X} 


1/2 


for  p v V,  by  (1.1),  (.1.1)  and  (1.7).  Result  (1.6)  remains  true  for  > (p-l)  v 1, 

and  If  (0v(p-l))<a  <l,  then  Rtp.u^.ail^)  * 0. 

Clearly,  R ' 1 if  1'  ' 2 and  R < 1 if  p »1 . Set 


p*  - (a0-l  + 2aP0)  («0-ltaBQ)  (1  < p1 


2) 


Umple  ulqobra  qives  R *=  1 if  p * p*,  and  R s 1 (K  - 1)  if  p ' p*  (p  ■ p*)  . It 


1 v p 


and  u - 1 , set 


1 + aB  (2-P)  (p-l)”'  . and 


(aB0)*  - l''0~  " (P-l)  < 2-p) 


-1 


Simple  ilqebra  givmi  R - 1 if  n * a*  , equivalently  aB^  - (aB^)  * , and  R s l (R  ' 1)  it 
a i*  (.<0  • n»)  or,  equivalently,  aB()  <■  (aB^l*  (a|l()  •>  (aB())M.  Finally,  i«  i easy  to  see 
that  R takes  on  all  values  In  |0,«)  . 

Til  summary,  for  qiven  values  of  p,  cft,  c^,  and  B(1,  Hi  and  11.'  can  l«  • 'mpnicd 

readily,  elfhet  i ■ inn  exact  values  of  V.  and  V,  or  an  asymptotic  approximation 


-u- 


I . . ^ 


E(C(t)).  The  exact  value 


Section  4.  An  Asymptotic  Comparison  of i with  B1  and  B2 

Analogous  to  the  earlier  definitions  of  and  V,,  let  VQ  * 

of  V()  must  be  computed  numerically,  but  its  limiting  value  is  given  in  Theorem  2.3.  The  class 
of  sequential  procedures  includes  B1  and  B2i  hence,  1 im  V()V.*  < 1,  for  i » 1,2,  with  the 
limit  taken  as  sampling  costa  tend  to  tern.  Thus,  baaed  on  the  criterion  of  expected  total  cost, 
neither  HI  nor  B2  are  ever  better  than  t in  the  limit.  However,  in  real  applications, 
other  considerations  t such  as  case  of  implementation)  are  important  and  an  experimenter  may 
prefer  HI  < H J ) if  ia  neai  unity.  In  this  section  the  limitinq  value  of  * , 

l - 1,2,  will  be  given  in  the  case  of  exactly  one  sampling  cost  positive.  In  the  case  of  both 
sampling  costs  positive  the  limit  of  v must  be  computed  numerically]  hence  it  is  not  a con- 
venient case  to  determine  general  patterns. 

It  ia  not  difficult  to  see  that  the  conclusions  of  Theorem  2.2  and  2.3  remain  true  if  the 
hypotheses  are  weakened  to  up  ■*  0 v (p-i).  In  fact,  if  * 0,  then  the  hypotheses  can  be 

weakened  to  > 0 v (p-2) . Similarly,  results  given  on  the  limiting  form  of  V,  can  be  ex- 
tended to  ci()  > j v (p-2),  where  j - 1 if  c > 0 and  j * 0 if  c^  « 0.  As  mentioned  in 
Section  3,  results  given  on  the  limitinq  form  of  require  only  ot^  > 0 v (p-1)  . 

For  -»>  v r < >”  and  b > 0 v (r-1)  , define 

H(b,r)»r(b+  (l-r)/2)i'(b)'l/2r<b+l-r)'1/2  . 


For  i “ 1,2,  set 


Q (T,i)  « lim  VV 

P C *0  ° 1 

A 

c -0 
T 


-1 


and  define  ( A ,1 ) in  the  analogous  way.  It  follows  easily  from  Theorem  2.2  and  2.1,  formu- 
las (1.7),  (3.1)  and  (3.3)  and  the  remarks  above  t.hat 


0|tiA,  l) 

- H(n0,p) 

’ ao 

> 

0 

V 

(p-D 

0p(A,2) 

- H(o0-l,p-2) 

no 

N 

l 

V 

(P-1) 

- 0 

1 i ®o 

S 

0 

V 

(p-i) 

- H(u0+l,p+l) 

n0 

0 

V 

p-i 

Q (T.2) 

- H(q0,p-1) 

a0 

> 

0 

V 

p-2 
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14. 1) 


Some  Examples  of  (4.1) 


for  simplicity,  only  compare  r to  the  "better"  of  B1  and  B2  (i.e.  B1  if  p ' 1. 


B2 

if 

p > 2,  see  Section  3) 

KJL 

_0: 

In  this  case 

Q0(T,1)  -1  , aQ  > 0,  and 

Qq ( A . 1 ) - r<a0  ♦ i)r(a0)_X/2r(a041)‘1/2,  a0  > 0 . 

Note  lim  Q(A,l)  - 0 and  lim  Q(A,1)  - 1. 

V° 

p - 1 1 In  this  case, 

O^A.l)  - 1 , 

* 

O^T.l)  - r<o0  4 i-)r(a0)"1/2r(a04l)"i/2,  c0  > 0 , 

and  the  remarks  of  the  previous  case  apply. 
p ■ 2 ; In  this  case, 

02(T,2)  - 1 , 

Q2(A,2)  - r(aQ  - j)r(a0-l)'l/2r(a0)'1/2,  aQ  > 1 , 

and  the  remark  of  case  p * 0 apply  with  the  modification  4 l instead  of  o»^  4 0. 
p *• 3s  In  this  case, 

Q3(A,2>  * 1 

Q,(T,2)  - P(a0  - i)r<rt0-l)'1/2r(n0)'l/2,  nQ  > l , 

and  the  previous  remark  applies. 
p < 0;  In  this  case, 

<yr,n  - r<n04i-p/2)r<o04i)'1/2r<«04i-p)'1/2,  and 
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